PARTICLE PRODUCTION IN HIGH-ENERGY HEAVY-ION 

COLLISIONS 



XIN-NIAN WANG 
Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, 

California, 94720 

Department of Physics, Shandong University, Jinan, P.R.China, 250100 

Particle production mechanisms in high-energy heavy-ion collisions are reviewed 
in connection with recent experimental data from RHIC. Implications on mini-jet 
production, parton saturation and jet quenching are discussed. 



1 Introduction 

Quark-gluon plasma (QGP) is expected to form in high-energy heavy-ion col- 
lisions. Though some proposed signals can provide more direct measurements 
of the initial parton density, many must reply on information inferred indi- 
rectly from the measurement of final hadron spectra. Therefore, both global 
observables such as the rapidity density of hadron multiplicity and high px 
hadron spectra can provide important links of a puzzle that can eventually 
lead one to a more complete picture of the dynamics of heavy-ion collisions 
and formation of QGP. 

Mini-jet production in a two-component model hasJoirg been proposed 
to explain, the energy dependence of total cross section Em and particle, pro- 
duction Era in high-energy hadron collisions. It has also been proposed Uu and 
incorporated in the HIJING model Lffl to describe initial parton production 
in high-energy heavy-ion collisions. In this model, parton production from 
semi- hard processes becomes quite significant, giving rise to the high initial 
parton density in high-energy heavy-ion collisions. However, there are still 
quite large uncertainties due to the lack of knowledge of the initial parton 
distributions in nuclei and the final state interaction of produced partons. 
Therefore, the study of energy and centrality dependence of central rapid- 
ity density □ can provide important constraints on models of initial entropy 
production and shed lights on the initial parton dffiphiotkms in nuclei. For 
example, the available RHIC experimental data li3'E3E3llj'L3 can already rule 
out the simple two-component model without nuclear modification of the par- 
ton distributions in nuclei □. However, further studies are needed to constrain 
the unknown nuclear shadowing of the gluon distribution in nuclei and to 
further distinguish the conventional parton production mechanism from other 



novel physics such as parton saturation 



s. 



Among the produced particle, a few have very high transverse momen- 
tum. They are related to high-energy jet production which is a very useful 
tool to study the properties of dense matter formed in high-energy heavy-ion 
collisions. Because large pr partons are produced very early in heavy-ion col- 
lisions and their production rates can be calibrated in pp and pA collisions at 
the same energy, they are ideal probes of the dense matter that is formed in 
the same reaction. What probes the dense medium is the scattering induced 
energy loss suffered by an energetic parton as it propagates through the mat- 
ter. The parton energy loss is directly related to the parton density of the 
medium. These are the main topics of this talk. 



2 Two-component model 

In the two-component model, the soft and hard processes are separated by a 
cut-off scale p . While the cross section of soft interaction cr soft is considered 
nonperturbative and thus noncalculable, the jet production cross section <7j et 
is assumed to be given by perturbative QCD (pQCD) for transverse momen- 
tum transfer px > po- The two parameters, a so ft and po, are determined 
phenomenologically by fitting the experimj^aLdata of total p + p(p) cross 
sections within the two-component model El'HaQLl'Q. The cut-off scale po, sep- 
arating nonperturbative and pQCD components, could in principle depend on 
both energy and nuclear size. Using the Gluck-Reya-Vogt (GRV) parameter- 
ization til of partem distributions and following the same procedure as in the 
original HIJING 0, one finds that the cut-off scale must have a weak energy 
dependence, 

p (Vs) = 3.91 - 3.34 log(logVi) + 0.98 log 2 (log y/s) + 0.23 log 3 (log y/s), (1) 

in order to fit the experimental data, because of the rapid increase of gluon 
distribution at small x. Shown in Fig. [I] is the calculated central rapidity 
density, 

dN ch Cj e t(s) 

— — - (n) s + ( n ) h ^-_w (2) 
drj Cin(s) 

for p + p(p) collisions as a function of energy y/s, where (n) s = 1.6 and (n)h = 
2.2 represent particle production from soft interaction and jet hadronization, 
respectively. The jet cross section is the lowest order pQCD result with a 
if- factor of 2 and GRV parton distributions. 

To extrapolate the two-component model to nuclear collisions, one as- 
sumes that multiple mini-jet production is incoherent and thus is proportional 
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Figure 1. Charged parfcwle-irapidity density per participating nucleon pair versus the cm. 
energy. The RHIC data f£Hlj (filled circle and ujr^rtarrale) for the 6% most central Au+Au 
are campared to pp and pp data (open symbols) ESEjEIj and the NA49 Pb + Pb(central 5%) 
data Ej (filled square). The two-component mini-jet model with and without shadowing 
is also shown. The shaded area for central Au + Au collisions corresponds to the range of 
gluon shadowing parameter s g = 0.24—0.28 [Eq. (kl)]. 



to the number of binary collisions Aginary- The soft interaction is however co- 
herent and proportional to the number of participant nucleons -/Vp art accord- 
ing to the wounded nucleon model E3. One important nuclear effect we have 
to consider in the two-component model is the nuclear shadowing of parton 
distributions in nuclei at small x. Such a nuclear shadowing effect in jet pro- 
duction can be taken into account by assuming modified parton distributions 
in nuclei, 

tf(x, Q 2 ) = AR*{x, Q 2 )f»(x, Q 2 ). (3) 

Using the experimental data from DIS off nuclear targets and unmodified 
DGLAP evolution equations, one can parameterize R^(x,Q 2 ) for different 
partons and nuclei. Recent new data however indicate that the simple, pa- 
rameterization for nuclear shadowing used in HIJING is too strong t3 for 
heavy nuclei. Instead one has to use the following new parameterization, 

(x) = 1.0 + 1.19 log 1/6 A (a; 3 - 1.2x 2 + 0.21a;) 



- s a (A 1 / 3 - 1)°' 6 (1 - exp(-a; 2 /0.01) (4) 



with s q = 0.1 for all quark distributions and s g = 0.24-0.28 for gluon. Assum- 
ing no final state effects on multiplicity from jet hadronization, the rapidity 
density of hadron multiplicity in heavy-ion collisions as shown in Fig. [l] is, 

^ = \(N pmt )(n) s + (n) h (N hinaiy )^^, (5) 

0,7] Z (T[ n 

where <r^(s) is the averaged inclusive jet cross section per NN in AA colli- 
sions including parton shadowing effect. The average number of participant 
nucleons and number of binary collisions for given impact-parameters can be 
estimated using HIJING Monte Carlo simulation, ft is clear that one has to 
consider parton shadowing in jet production in heavy-ion collisions. 



3 Parton Saturation 



Similar results of multiplicity in heavy- ion collisions are also predicte d by other 
models EjO, in particular the initial-state partorLsaluration model E£lil3. It is 
based on the nonlinear Yang- Mills field dynamics Bt3 assuming that nonlinear 
gluon interaction below a saturation scale Q 2 S ~ a s xGa(x,QI)/ttR a leads to a 
classical behavior of the gluonic field inside a large nucleus, where Ga(x, Q 2 s ) 
is the gluon distribution at x = 2Q s /*/s. Assuming particle production in 
high-energy heavy-ion collisions is dominated by gluon production from the 
classical gluon field, one has a simple form Il3 for the charged hadron rapidity 
density at 77 = 0, 



dN, 



ch 



(ATpart) drj 



so 



A/2 



log 



A 2 



— log f — 

2 S Vso 



(6) 



with c w 0.82 111. This is shown in Fig. |^ as solid lines as compared to the two- 
component model (shaded area). Here, Aqcd = 0.2 GeV, A = 0.25 and the 
centrality jdependence of the saturation scale Qq s at ^/so = 130 GeV is taken 
from Ref. cB. The two-component model has an impact-parameter dependent 
parton shadowing. 

Comparing the two model results in Fig. |2|, one notices that the saturation 
and two-component model agree with each other in most regions of centrality 
except very peripheral and very central collisions. In central collisions, re- 
sults of saturation model tend to be flatter than the two-component model. 
In this region, there are still strong fluctuations in parton production in the 
two-component model through the fluctuation of iVbi nar y while N pax t is lim- 
ited by its maximum value of 2 A. That is why dN c h/ drj/ (N par t) continues to 




Figure 2. The charged hadron central rapidity density per participant nucleon pair as a 
function of the averagedrf*umber of participants from the two-component imodel (shaded 
lines), two-parameter [fit pi -(«ieft|-dashed lines) and parton saturation model E£| as compared 
to experimental data H3'E3li£pE3. 



increase with (-/V part ) in the central region. Such a fluctuation is not currently 
taken into account in the saturation model calculation. More accurate mea- 
surements with small errors (less than 5%) will help to distinguish these two 
different behaviors. For peripheral collisions, saturation model results fall off 
more rapidly than the mini-jet results. However, the experimental errors are 
very big in this region because of large uncertainties related to the determi- 
nation of the number of participants. Therefore, it will be very useful to have 
light-ion collisions at the same energy to map out the nuclear dependence of 
the hadron multiplicity in this region. An alternative is to study the ratios of 
hadron multiplicity of heavy-ion collisions at two different energies as a func- 
tion of centrality c3. In this case, the errors associated with the determination 
of centrality will mostly cancel. 

It is interesting to point out that in the saturation model that assumes 
a particle production mechanism dominated by coherent mini-jet production 
below the saturation scale Q s , the value of Q s determined in Ref.li3't£l is 
much smaller than the cut-off po in the two-component model constrained by 
the p + p(p) data. As demonstrated in this paper, the number of mini-jet 
production below such scale is still very large and should contribute to the 



final hadron multiplicity. 



4 Parton energy loss 

Among the bulk of produced hadrons in heavy-ion collisions, there are also 
high pt hadrons produced in events with hard processes. In principle, jet 
quenching effect could also lead to increased total hadron multiplicity c3_duc 
to the soft gluons from the bremsstrahlung. However, a recent study c3 of 
parton energy loss in a thermal environment found that the effective energy 
loss is significantly reduced for less energetic partons due to detailed balance 
by thermal absorption. Thus, only large energy jets lose significant energy via 
gluon bremsstrahlung. Since the production rates of these large energy jets 
are very small at the RHIC energy, their contributions to the total hadron 
multiplicity via jet quenching should also be small. Similarly we also assume 
that parton thermalization during the early stage contributes little to the final 
hadron multiplicity. 

Theoretical studies of the parton energy loss in hot medium date back 
to the first attempt by Bjorken E3 to calculate elastic energy loss of a parton 
via elastic scattering in the hot medium. A simple estimate can be given 
by the thermal averaged energy transfer v e \ w q^_/2ui of the jet parton to a 
thermal parton with energy ui, q± being the transverse momentum transfer of 
the elastic scattering. The resultant elastic energy loss 13 



dx 2 V V 

is sensitive to the temperature of the thermal medium but is in general small 
compared to radiative energy loss. Here, fi is the Debye screening mass and 
Ci is the Casimir of the propagating parton in its fundamental presentation. 
jThe elastic energy loss can also be calculated within finite temperature QCD 
E3 with a similar result, but with a more careful and consistent treatment of 
screening effect. 

Though there had been, estimates of the radiative parton energy loss using 
the uncertainty principle till, a first theoretical study of QCD radiative parton 
energy loss incorporating Landau-Pomeranchuk-Migdal interference effect 123 
is by Gyulassy and myself c3 where multiple parton scattering is modeled by a 
screened Coulomb potential model. Baier et al. (BDMPS) Lj later considered 
the effect of gluon rescattering which turned out to be very important for 
gluon radiation induced by multiple scattering in a dense medium. These two 
studies have ushered in many recexdj works on the subject, including a path 
integral approach to the problem Ej and opacity expansion framework LjO 



which is more suitable for multiple parton scattering in a thin plasma. The 
radiative parton energy loss to the leading, order of the opacity n = L/X in 
the thin plasma of size L is estimated as ESEa 

dE iad a sf i 2 L (2E\ 

where A is the gluon's mean-free-path in the medium. The unique L depen- 
dence of the parton energy loss is a consequence of the non-Abelian LME 
interference effect in a QCD medium. It is also shown in a recent study E3 
that thermal absorption and stimulated emission in a thermal environment 
can be neglected for high energy partons (E 3> n) while they are important 
for intermediate energy partons. 

Using this latest result one can estimate the total energy loss for a parton 
with initial energy E = 40 GeV to be about AE « 10 GeV after it propagates 
a distance of L = 6 fm in a medium with \i = 0.5 GeV and A = 1 fm. For 
an expanding system—the total energy loss is reduced by a factor of 2tq/L 
from the static value Assuming that most of this energy loss is carried 

by gluons outside the jet cone 113, measuring the energy loss would require 
the experimental resolution 5E to be much smaller than the total energy loss 
AE. With the measured totaljiiultiplicity density dN/drj « 900 t2l and energy 
density dEx/dr] « 500 GeV eil in central Au + Au collisions at \fs = 130 
GeV, one can estimate that the average total background energy within the 
jet cone (Sj] = 1 and 8<fr = 1) is about HEt ~ 80 GeV with a fluctuation of 
SEt ~ 10 GeV. It is therefore very difficult, if not impossible, to determine the 
energy of a jet on a event-by-event base c3. Since high px hadrons in hadron 
and nuclear collisions come from fragmentation of high Et jets, energy loss 
naturally leads to suppression of high pr hadron spectra. Miklos Gyulassy 
and I then proposed E3 that one has to reply on measuring the suppression 
of high pt hadrons to study parton energy loss in heavy-ion collisions. Since 
inclusive hadron spectra is a convolution of jet production cross section and 
the jet fragmentation function in pQCD, the suppression of inclusive high 
Pt hadron spectra is a direct consequence of the medium modification of the 
jet fragmentation function induced by parton energy loss. Assuming that jet 
fragmentation function is the same for the final leading parton with a reduced 
energy, the modified fragmentation function can be assumed as l3 

- T-Ke!e d {i^eJe) • ^ 



Therefore, in this effective model the measured modification of fragmentation 
function can be directly related to the parton energy loss. 



5 Modified fragmentation functions 

Since a jet parton is always produced via a hard process involving a large 
momentum scale, it should also have final state radiation with and without 
rescattering leading to the DGLAP evolution equation of fragmentation func- 
tions. Such final state radiation effectively acts as a self-quenching mechanism 
softening the leading parton momentum distribution. This process is quite 
similar to the, induced gluon radiation and the two should have strong inter- 
ference effect ESEZJ. It is therefore natural to study jet quenching and modified 
fragmentation function in the framework of modified DGLAP evolution equa- 
tions in a medium CJ. 

The simplest case for jet quenching is deeply inelastic scattering of an 
electron off a nucleus target where the virtual photon knocks one quark out of 
a nucleon inside a nucleus. The quark then will have to propagate through the 
rest of the nucleus and possibly scatter again with other nucleons with induced 
gluon radiation. The induced gluon radiation reduces the quark's energy 
before it fragments into hadrons with a modified fragmentation function. One 
can study the nuclear modification of the fragmentation function by comparing 
it with the same measurement in DIS with a nucleon target. 

In an infinite momentum frame, where the photon carries momentum q = 
[—x bP + , q~ , Ox] arid the momentum of the target per nucleon is p = [p + , 0, 0_J 
with the Bjorken variable defined as xb — Q 2 /^q p + , one can calculate the 
modified fragmentation function. Because of the LPM interference effect, the 
formation time of the gluon radiation due to the LPM interference requires the 
radiated gluon to have a minimum transverse momentum ~ Q 2 /MRa ~ 
Q 2 /A 1 / 3 . The nuclear corrections to the fragmentation function due to double 
parton scattering will then be in the order of otsA 1 / 3 /^ ~ a s A 2 l 3 /Q 2 , which 
depends quadratically on the nuclear size. For large values of A and Q 2 , these 
corrections are leading and yet the requirement i\ <C Q 2 for the logarithmic 
approximation in deriving the modified fragmentation function is still valid. 

Shown in Fig. || are the calculated nuclear modification factor of the 
fragmentation function for 14 N and 84 ifr targets as compared to the recent 
HERMES datac3. The predicted shape of the z dependence and the quadratic 
nuclear size dependence agrees well with the experimental data. The energy 
dependence of the suppression also has excellent agreement with our predic- 
tion. What is amazing is the clear quadratic A 2 ! 3 nuclear size dependence of 
the suppression which is a true QCD non-Abelian effect. In fitting the data 
of the overall suppression for 14 iV target we obtain the only parameter in our 
calculation, Ca 2 = 0.00065 GeV 2 . 

One can also calculate theoretically the average energy loss by the quark, 




which is the energy carried away by the radiated gluons, 

AE = v(Az g ) « C(C A a 2 /N c )3ln(l/2x B ). (10) 

With the value of a 2 C, one gets the quark energy loss dE/dx « 0.3 GeV/fm 
for the Kr target. 

6 Jet quenching in heavy-ion collisions 

In high-energy heavy-ion collisions, the jet production rate is not affected by 
the formation of dense matter and the final state multiple scattering. One can 
assume that the high px hadron spectra can then be given by the convolution 
of the jet production cross section and the medium modified jet fragmentation 
function D h / c {z c ,Q 2 ). 

In principle, one should use the modified fragmentation function evaluated 
according to the pQCD calculation for a dense medium. However, before that 
can be done in a practical manner, we have used the effective approach in 
Eq. (||) by rescaling the fractional momentum by 1 — Az to take into account 
of the parton energy loss. To verify whether such an effective approach is 
adequate, we compare the two modified fragmentation functions in Fig. ^. 
We found that the effective model (dashed lines) can reproduce the pQCD 
result (solid lines) very well, but only when Az is set to be Az w 0.6(z g ). 
Therefore the actual averaged parton energy loss should be about 1.6 times 
of that used in the effective modified fragmentation function. This difference 




is caused by the absorptive processes or unitarity correction effect in the full 
pQCD calculation. 

Unlike in D1S nuclear scattering, the dense medium in high-energy heavy- 
ion collisions is not static. It has to go through rapid expansion which should 
also affect the effective total parton energy loss. The total energy loss ex- 
tracted from experiments should be a quantity that is averaged over the whole 
evolution history of the expanding system. It is therefore useful to convert the 
averaged quantity to an energy loss in a static system that has the same par- 
ton density as the expanding system at its initial stage. If the averaged total 
parton energy loss in aJongitudinally expanding system with a transverse size 
R is AEid, one finds t£l that the corresponding parton energy loss in a static 
system with the same initial parton density would be AE = AEi < \(R/2tq). 
Here tq is the initial formation time of thc_dense medium. 

Comparing the recent PHENIX data c3 with pQCD parton model calcu- 
lation as shown in Fig. [|, one can extract a value of dE/dx=0.25 GeV/frn that 
one needs to use in the effective modified fragmentation function in fitting the 
data. Taking into account the unitarity correction effect and the expansion, 
this corresponds to an effective energy loss dE/dx = 1.6 X 0.25R/2tq in a static 
system with a density similar to the initial stage of the expanding system at 
To. With R ~ 6 fm and To ~ 1/po = 0.2 fm, this would give dE/dx « 12 
GeV/fm, which is about 40 times of that in a cold nuclear matter as extracted 
from the DIS processes. Since the parton energy loss is directly proportional 
to gluon density, this implies that the gluon density in the initial stage of 
Au + Au collisions is about 40 times higher than that inside a cold nucleus. 
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Figure 5. pQCD parton model calculation of the charged hadron andr-njmi spectra in pp 
and central Au + Au collisions compared with the experimental data HScJ. The effective 
modified fragmentation function is used in the calculation. 



7 Conclusions 

Recent RHIC data require a strong shadowing of gluon distribution in nuclei 
within a two-component model of particle production in heavy-ion collisions. 
Using this strong gluon shadowing with an assumed impact-parameter depen- 
dence, the predicted centrality dependence of the hadron multiplicity agrees 
well with the recent. RHIC results. The results are compared with the parton 
saturation model lia-tll In order to differentiate the two models one needs 
more accurate experimental data in both the most central and peripheral re- 
gions of centrality or study the centrality dependence of the ratios at different 
colliding energies. 

It is argued that the jet quenching will not contribute much to the total 
hadron multiplicity due to detailed balance in a thermal environment, even 
though they significantly suppress the high px spectra. The medium modifica- 
tion of the jet fragmentation functions due to gluon radiation induced by the 
multiple parton scattering can be calculated in pQCD. The predictions of the 
shape, energy dependence, and the quadratic nuclear size A 2 / 3, dependence 
of the modification agree well with the recent HERMES data. The resultant 
parton energy loss in the cold nuclear medium is estimated to be about 0.3 



GeV/fm inside Kr nuclei. Comparing to the QCD result of the modifica- 
tion of fragmentation function, the actual averaged energy loss is found to 
be about 1.6 times that of the effective energy loss used in a earlier effective 
model for the same modification. Considering the effect of expansion, the 
recent PHENIX data imply a medium induced energy loss in central Au + Au 
collisions equivalent to 12 GeV/fm in a static medium with the same gluon 
density as in the initial stage of the collision. 
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